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Abstract. We present the analysis of a large sample of lines of sight with extinction curves 
covering wavelength range from near-infrared (NIR) to ultraviolet (UV). We derive total to 
selective extinction ratios based on the Cardelli, Clayton & Mathis (1989, CCM) law, which 
is typically used to fit the extinction data both for diffuse and dense interstellar medium. We 
conclude that the CCM law is able to fit most of the extinction curves in our sample. We divide 
the remaining lines of sight with peculiar extinction into two groups according to two main 
behaviors: a) the optical/IR or/and UV wavelength region cannot be reproduced by the CCM 
formula; b) the optical/NIR and UV extinction data are best fit by the CCM law with different 
values of Rv- We present examples of such curves. The study of both types of peculiar cases 
can help us to learn about the physical processes that affect dust in the interstellar medium, 
e.g., formation of mantles on the surface of grains, evaporation, growing or shattering. 



1. Introduction 

Interstellar grains affect starlight which passes through them by absorbing and scattering 
photons. These two physical processes produce the interstellar extinction which depends on 
the properties of dust grains, e.g., size distribution and composition. The description of the 
average extinction curve of our Galaxy as a function of the wavelength from the infrared to 
ultraviolet can be found in Savage pQ. Extinction curve shows some evident features: it rises in 
the infrared, it shows a slight knee in the optical, it is characterized by a bump at 2175^4, and 
it rises in the far- ultraviolet. These features are common between different environments. The 
properties of interstellar grains are different in diffuse and dense interstellar medium and thus 
also the extinction changes. CCM |2 a found an average extinction law valid over the wavelength 
range 0.125 fim < A < 3.5/j.m, which is applicable to both diffuse and dense regions of the 
interstellar medium. This extinction law depends on only one parameter Ry = Ay/E{B — V). 
The Ry parameter ranges from about 2.0 to about 5.5 (with a typical value of 3.1) when one goes 
from diffuse to dense interstellar medium and thus Ry characterizes the region that produces 
the extinction. 

If one knows the value of Ry along a particular line of sight, one can obtain the extinction 
curve from the infrared to ultraviolet using the CCM law (see Figure ^): 

4^ = a(x)+b(x)-R v \ (1) 




where x = 1/A, and a(x) and b(x) are the wavelength-dependent coefficients. 

There are different ways to obtain Ry using the NIR or UV extinction data. Wegner [2] 
computed Ry values for the sample of 597 OB stars with known near- infrared magnitudes. He 
assumed that, in the infrared spectral region, the normalized extinction curve is proportional to 
A -3 or A~ 4 . Extrapolating the IR interstellar extinction curve to 1/A = he derived Ry as: 



Ry 



E{\ - V) 
E(B - V) 



(2) 



Gnacihski & Sikorski 4 applied the x 2 minimization method to compute the Ry values for a 
sample of ultraviolet (UV) extinction data using the linear relation Geminale & Popowski 
jS] extended the analysis from 0j by using non-equal weights derived from observational errors 
to determine Ay and Ry values toward a sample of stars with known ultraviolet color excesses. 

In this paper we use both optical/infrared (URIJHKLM) and UV extinction data to obtain 
Ry values for a sample of 436 lines of sight. We arrive at two main conclusions: (i) there are lines 
of sight with extinction in the optical/IR or/and UV which generically don't follow the CCM 
law, and (ii) there are lines of sight which show an extinction curve that cannot be reproduced 
with a single Ry value in the whole wavelength range. 

2. Theoretical Basis 

We normalize the extinction in a standard way: 



(3) 



The absolute extinction may be deduced from the relative extinction by using a total-to-selective 
extinction ratio Ry. 

Rv = ^TTT- (4) 



Then: 
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For each individual band, equation Q and Q can be combined to derive an Ry value. More 
generally, the \ 2 minimization can be used to obtain the Ry value that provides the best CCM 
fit to all observed extinction data. Gnacihski & Sikorski 4 suggested the following \ 2 > which 
can be minimized to derive Ry. 
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where a(xi) and b{xj) are the coefficients of the CCM law. Following our previous work we 
use an improved weighted formula to find Ry. We minimize the following \ 2 '- 
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where u>i = 1/crf are the weights associated with each band. 
Minimizing equation (JJJ) with respect to Ry, we find: 
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where here Oi values are taken from Wegner Q and they were computed according to: 
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We estimate the error in Ry from: 
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3. Data 

We use a sample of 436 lines of sight with optical/IR and UV extinction data reported by 
Wegner (HJ. In Wegner's computation the ultraviolet photometry is taken from Wesselius et al. 
[7] and based on Astronomical Netherlands Satellite (ANS); infrared magnitudes in J, H, K, L, M 
passbands are originate mostly from the catalog of Gezari, Schmitz and Mead !8 and Gezari et 
al. [S]. The R and I magnitudes with accuracy of ±0.01 mag are taken from Johnson ^Uj and 
Fernie The spectral classification and UBV data which have the accuracy of ±0.01 mag 
come from the SIMBAD database. 

The effective wavelengths of the optical/IR bands are: Xjj = 0.36/um, Xr = 0.71/um, 
Aj = 0.97/im, Aj = 1.25/xm, Xh = 1.65/xm, Xk = 2.2/xm, Al = 3.5/xm, Xm = 4.8//m; whereas 
the effective wavelengths of the UV bands are: 0.1549, 0.1799, 0.2200, 0.2493, and 0.3294 fim. 

From the original sample, we exclude 20 lines of sight because of their negative value of 
Ry. Negative Ry values are typically the result of noisy data for the lines of sight with small 
E(B-V). 
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Figure 2. Examples of extinction curves with x 2 /dof > 2.0. The long dashed line represents 
the CCM curve which is the best fit to the IR data and the point-short dashed line represents 
the CCM curve which best follows the UV data. 



4. Results 

4.1. x 2 test 

A useful method to test if all extinction data points are well fitted by the CCM law is to compute 
X 2 based on equation normalized to the number of degrees of freedom, which is equal to 
the number of observed points minus the number of fitted parameters (in our case the only 
parameter is Ry). We note that our average x 2 /dof is not equal to the expected value of 1; 
since the errors are taken from Wegner (2002) and we do not want to modify them, we do not 
renormalize our errors requesting < x 2 /dof >= 1. Therefore it is safer to treat our x 2 /dof 
as a measure of a relative rather than absolute quality of the CCM fit in the optical/IR and 
UV wavelength ranges. We select outliers considering the tail of our x 2 /dof distribution. We 
assume that the values x 2 /dof > 2.0 are indicative of the lines of sight with extinction curves 
not well fitted by the CCM law. We find that 25% of the lines of sight of our sample shows this 
disagreement with the CCM law. Figure |2 displays two extinction curves representative of this 
group. 

4-2. Test for the universality of CCM law 

We further analyze the lines of sight with x 2 /dof < 2.0 for which we expect the CCM law to 
fit well all observed extinction data from optical/IR to UV. The usual assumption is that the 
knowledge of the Ry value obtained from the IR part of the extinction curve may be used to 
obtain the entire extinction curve by using the CCM law. We critically test this assumption 
using two sets of Ry values for each line of sight: the infrared Ry ( lR Ry) and ultraviolet Ry 




Figure 3. The comparison between the Ry values obtained from optical/IR and UV extinction 
data. The points with error bars are those for which \5\ > 2.0. 



( Ry). We compute the following statistic: 

^Ry - lR R, 



a 2 [^R v ]+a 2 [ 1R Ry] 



(11) 



Specifically, when \5\ > 2.0 we classify the line of sight as anomalous in the sense that the CCM 
law is not able to reproduce the whole extinction curve with a single value of Ry. 

Figure shows the comparison between the Ry values obtained from optical/IR and UV 
extinction data. There is a large scatter around the 1-to-l relationship; however, only the points 
with the error bars shown are those which deviate from this relation significantly (\5\ > 2.0). 
These 25 lines of sight need two different values of Ry to match the optical/IR and UV part of 
the extinction curve properly. Therefore, in these cases we cannot universally use the CCM law 
to represent the entire extinction curve. Six cases out of 25 have the uv iZy higher than the 1R Ry 
one, and the other 19 have the uv i?y lower than the one. Figure 0] shows two examples of 
such anomalous extinction curves. 



5. Conclusion 

We use a x 2 minimization method to compute the Ry values for a sample of 436 lines of sight. 
We exclude 20 cases because of their negative Ry and analyze the final sample of 416 lines of 
sight. We derive our Ry values assuming CCM law. This law aims at reproducing the entire 
extinction curve by using only one parameter: Ry. We analyze the goodness of the CCM fit 
for all lines of sight using the x 2 /dof statistic and test the universality of the CCM law by 
computing Ry values separately for the optical/IR and UV part of the extinction curve. We 
find that for 69% of our original sample, the CCM law is able to fit well both the optical/IR 
and UV data. We divide the remaining 31% of cases into two groups, according to two main 
peculiarities: a) the optical/IR or/and UV extinction data points cannot be fitted by the CCM 
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Figure 4. Examples of extinction curves for which two different Ry values are necessary to 
obtain a good fit to the entire extinction curve. 



law (25% of the entire sample); b) Ry values are significantly different for the two spectral 
regions: optical/IR and UV (6% of the entire sample). Unless caused by faulty data, peculiar 
extinction curves result from unusual properties of dust grains. Therefore, theoretical modeling 
of these extinction curves (e.g., Mishchenko ^2j; Saija et al. |13j : Weingartner & Draine [Tl 
may help us to understand the processes which modify the properties of interstellar grains. 
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